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ABSTRACT 

In this work, we examine the effects of mergers on the hydrostatic mass estimate of galaxy clusters 
using high-resolution Eulerian cosmological simulations. We utilize merger trees to isolate the last 
merger for each cluster in our sample and follow the time evolution of the hydrostatic mass bias 
as the systems relax. We find that during a merger, a shock propagates outward from the parent 
cluster, resulting in an overestimate in the hydrostatic mass bias. After the merger, as a cluster 
relaxes, the bias in hydrostatic mass estimate decreases but remains at a level of -5-10% with 15-20% 
scatter within r^. We also investigate the post-merger evolution of the pressure support from bulk 
motions, a dominant cause of this residual mass bias. At rsoo, the contribution from random motions 
peaks at 30% of the total pressure during the merger and quickly decays to ~ 10-1 5% as a cluster 
relaxes. Additionally, we use a measure of the random motion pressure to correct the hydrostatic 
mass estimate. We discover that 4 Gyr after mergers, the direct effects of the merger event on the 
hydrostatic mass bias have become negligible. Thereafter, the mass bias is primarily due to residual 
bulk motions in the gas which are not accounted for in the hydrostatic equilibrium equation. We 
present a hydrostatic mass bias correction method that can recover the unbiased cluster mass for 
relaxed clusters with 9% scatter at rsoo and 11% scatter in the outskirts, within r2oo- 
Subject headings: cosmology: theory galaxies: clusters: general 



1. INTRODUCTION 

Over the last decade, measurements of the abundance 
of galaxy clusters have begun to provid e important con- 
straints on cosmological parameters dVikhlinin et al.| 
2QQ9; |Mantz et al.|[2010l |Rozo et al.||2009| |Vanderlinde| 
et al. 2010; Sehgal et al.|2010| ). Clusters trace the highest 
peaks in the matter density field; their abundance as a 
function of mass and redshift is therefore extremely sen- 
sitive to the growth rate of density perturbations as well 
as the expansion rate of the Universe. 

The power of cluster number counts as a cosmological 
probe is principally limited by the accuracy with which 
their total mass can be measured. Cluster mass can be 
estimated using a wide variety of techniques, from mea- 
surements of the number or velocity dispersion of the 
cluster member galaxies, strong and weak gravitational 
lensing of background sources, and measurements of their 
X-ray emission or Sunyaev-Zel'dovich (SZ) signal. The 
two latter methods utilize the temperature and density 
of the hot plasma, or intra-cluster medium (ICM), that 
pervades cluster environments and dominates their bary- 
onic mass. The thermal properties of the ICM are very 
sensitive to the depth of the gravitational potential of the 
cluster host dark matter halo, and therefore to the total 
cluster mass. In particular, by measuring the radial den- 
sity and temperature profiles of the ICM from X-ray and 
SZ observations, the equation of hydrostatic equilibrium 
(HSE) can be solved to obtain the total matter profile. 

Application of the equation of HSE is based on the 
assumption that the ICM is both spherically symmetric 
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and in hydrostatic equilibrium. However, clusters are dy- 
namically young systems and obtain a significant fraction 
of their mass through mergers. Therefore the accuracy 
of this assumption must be questioned. Recent cosmo- 
logical simulations have suggested that even in relaxed 
clusters the HSE method may systematic ally underesti- 
mate the total cluster mass by ~ 10- 20% ( |Kay et aL 2004 
|Rasia et al.|20 06; Nagai et al. 2007b). It has been demon- 



strated that this bias is a consequence of neglecting the 
contribution of residual bulk motions in the ICM to the 



dynamical stability of relaxed clusters ([Kay et al.||"' 
Rasia et al.||2004[ IJeltema et al.||2008| IPitfaretti & 



larnini||2008| |Lau et al.||2009p. Typical 
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y, the HSE tec 
nique only accounts tor support from thermal pressure; 
the non-thermal contribution from bulk motions and tur- 
bulence is difficult to measure directly. The evolution of 
these bulk motions, their generation through mergers, 
and the potential effect on the corresponding evolution 
of the mass bias are all very poorly understood. 

Several recent works use cosmological simulations of 
clusters to investigate the effects of dynamical state on 
the deviation from hydrostatic equilibrium by examining 
a numb er of clusters at single redshifts. Nagai et al. 



(2007b| demonstrate that the mass bias increases rapidly 
with radius and is 20% or greater at all radii in unrelaxed 
systems. Utilizing a large sa mple of over 100 clusters, 
Piffaretti & Valdarnini (2008) show that the hydrostatic 



mass bias is larger tor the less relaxed, recently merged 
clusters, and also that the scatter in the bias increases 
with the degree of disturbance in the systems. 

In an attempt to incorpor ate observational uncertain- 
ties, [Meneghetti e t al.| ( |2010| ) investigate the effect of line 
of sight projections using mock lensing and mock X-ray 
observations and were able to recover the bias measured 
direct l y in s i mulations. In an e xtension of Meneghetti 
|et al.| fl2010| ), |Rasia et alT] ( |2012D also investigate the et- 
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fects of cluster morphology, environment, temperature 
inhomogeneity, and mass on the mass biases. Moreover, 
both works show agreement with an observational com- 
parison of the X-ray hydr ostatic mass estimate a nd weak 
lensing mass estimates by |Mahdavi et ah] ([2008 ) , namely 
that the HSE mass derived from Chandra X-ray mea- 
surements is biased low by ~ 20% at rsocQ compared 
to weak grav i tation al lensing measurements. Similarly, 



TABLE 1 
List of Simulated Clusters 



Zhang et al. (2010) investigate the effect of dynamical 
state on the bias between the mass measurement meth- 
ods and find a ~ 10% underestimate in the HSE mass 
at rsoo in undisturbed clusters but also find that in dis- 
turbed clusters the mass is overestimated by ~ 6% at the 
same radius. 

In this paper, we use high-resolution cosmological sim- 
ulations to investigate the time evolution of the hydro- 
static mass bias and non-thermal pressure support in 
the ICM as a function of time since the latest merger. 
While a number of recent studies have explored the ef- 
fects of mergers and dynamical state on the hydrostatic 



mass estimate (e.g. Kay et al. 
Nagai et al.||2007bl Lau et al. 



2004; |Rasia et al. 
2009; |Poole et al.' 



2006 
2006 



Piffaretti & Valdarnini||2008 ) , most only explored these 



effects at a single epoch or in idealized mergers. Our 
simulations and merger trees enable us to characterize 
the dynamical state of clusters robustly and quantita- 
tively, and more effectively follow the evolution of the hy- 
drostatic mass estimate during and after mergers in the 
context of realistic cosmological simulations. We investi- 
gate the source of the hydrostatic mass bias by following 
the post-merger evolution of the non-thermal pressure 
component of the ICM contributed primarily by random 
gas motion^] We find that during mergers, the bulk 
motions supply about 30% of the pressure support in 
clusters within rsoo- While the random motion pressure 
quickly decays, we also find that after 9 Gyr, these gas 
motions still provide about 10% of the pressure support. 
By taking the random motion pressure into account in 
a correction factor, we are able to recover the true mass 
of relaxed clusters, clusters which have not undergone a 
merger in more than 4 Gyr, to within a few percent. We 
conclude by presenting a potential method for estimating 
this correction factor using the thermal pressure profile. 

In Section [2] we describe our cosmological cluster sim- 
ulations. InSection 3} we present our analysis of the 
time evolution of the hydrostatic mass bias. In Section [4] 
we explore the connection between non-thermal pressure 
and present our model for correcting the mass bias. In 
Section [5| we summarize and discuss our findings. 

2. SIMULATIONS 

We analyze a sample of 16 simulated grou ps and clus- 



ters of gal axies presented pre viously in Na gai et al. 
( |2007b| ) and |Nagai et al.| ( [2007a| ), and we refer the reader 
to these papers tor more details. We briefly summarize 
the relevant parameters of the simulations here. 

The simulations are performed using the Adaptive 
Refinement Tree (ART) iV-body+gas-dynamics code 

4 Throughout this paper, we refer to total mass and ICM prop- 
erties within radii which correspond to fixed overdensities A rela- 
tive to the critical density at that redshift, such that Ma^a) = 
A(4/3)7Tr 3 A p c (z). 

5 Additional contributions from cosmic rays and magnetic field 
are not included. 
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1 Cluster labels correspond to those used in Nagai et al. (2007b). 

2 Time elapsed at z = since last major merger. 

3 At z = 0, CL107 is in the prelimary phases of a major merger 
with a halo that has not yet entered r2oo and therefore is classified 
as an unrelaxed cluster despite having tmerger = 5.98 Gyr. 



( |Kravtsov|[T999| |Kravtsov et~aT^|2QQ2l |Rudd et al.|2QQ8| > , 
which is an Eulerian code that uses adaptive refinement 
in space and time, and non- adaptive refinement in mass 
( Klypin et al.||2QQl| ) to achieve the dynamic range neces- 
sary to resolve the cores of halos formed in self-consistent 
cosmological simulations. The simulations assume a flat 
ACDM model: Q m = 1 - Q A = 0.3, Q h = 0.04286, 
h = 0.7 and os = 0.9, where the Hubble constant is 
defined as lOO/i km s _1 Mpc -1 , and as is the mass vari- 
ance within spheres of radius Sh~ r Mpc. The box size 
for CL101-CL107 is 120 h~ x Mpc comoving on a side and 
80/i _1 Mpc comoving for CL3-CL24. All simulations 
were run using a uniform 128 3 grid with 8 levels of mesh 
refinement, corresponding to peak spatial resolution of 
~ 7 h~ x kpc and 5 IrC 1 kpc for the two box sizes respec- 
tively. Only the regions ~ 3 — 10 /i — 1 Mpc surrounding 
each cluster center were adaptively refined. The dark 
matter (DM) particle mass in the region around each 
cluster was ra p ~ 9.1 x 10 8 h' 1 M for CL101-CL107 
and ra p ~ 2.7 x 10 8 h' 1 M for CL3-CL24, while other 
regions were simulated with a lower mass resolution. In 
order to isolate the effect of mergers on the ICM we use 
simulations th at treat on l y grav ity and non-radiative hy- 



drodynamics. |Lau et al. (20091) show that the exclusion 



of cooling and star formation does not have a strong ef- 
fect on the total contribution of gas motions to the HSE 
reconstruction of mass profiles outside of cluster cores. 

Each cluster has 30 complete outputs in the range 
z < 4, with the majority occurring equally spaced in ex- 
pansion factor after z = 1 with Aa = A(l+z) _1 = 0.025. 
The dark matter particles, used to identify cluster merg- 
ers, are output at a much finer spacing with ^300 outputs 
over the cluster history. This allows for a robust identi- 
fication of the cluster progenitor and merger timing. 

2.1. Halo Finder 

A vari ant of the bound de nsity maxima halo finding al- 
gorithm ( Kly pin et al.|1999| ) is used to identify halos and 
the subhalos within them. The details of the algorithm 
and parameters used in the halo finder can be found in 
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Kravtsov et al. (2004). The main steps of the algorithm 
are identification of local density peaks (potential halo 
centers) and analysis of the density distribution and ve- 
locities of the surrounding particles to test whether a 
given peak corresponds to a gravitationally bound clump. 
More specifically, we construct density, circular velocity, 
and velocity dispersion profiles around each center and 
iteratively remove unbound par ticles using the procedure 
outlined in Klypin et al. (1999). 



Following Lau et al. (2009), we treat the random gas 



2.2. Merger Trees 

To examine the evolution of cluster properties and the 
hydrostatic mass bias, we track the evolution of each 
cluster's most massive progenitor. We identify the main 
progenitor of each z — cluster by iteratively follow- 
ing the 10% most bound dark matter particles at each 
epoch. Mergers are located by identifying objects that 
cross r2oo of the parent cluster. We determine the mass 
ratio of the merger at an epoch when the two systems are 
sufficiently separated to avoid impacting their respective 
mass profiles, herein we use the epoch at which they ini- 
tially overlap at a radius of r2oo- We use the ratio of 
M 5 oo at this epoch to identify major mergers, herein de- 
fined as a mass ratio greater than 1:6. Merging events 
with the mass ratios smaller than this are defined as mi- 
nor mergers. While this r atio is somewhat lowe r than 
used in other studies (e.g. |Gottlober et al.||2001[ ), care- 
ful exploration of the merger ratio has shown that this 
threshold best separates mergers that have a strong im- 
pact on the measured ICM properties from those that 
do not. We examine the evolution of these properties as 
a function of time since the last major merger, t merger , 
defined as the epoch at which the center of the merging 
cluster first crosses r2oo of the parent cluster. 

2.3. Cluster Profiles 

For each simulation output we measure properties of 
the ICM under the assumption of spherical symm etry 
following the methods described in Lau et al. (2009). A 
brief description of the methods used is included here. 
We focus in this work on the estimate of cluster mass 
provided by the thermal hydrostatic equilibrium equa- 
tion, 

MHSE«r) = --A^, (1) 



Gp, 



gas 



dr 



where Mhse(< is the estimated mass within radius, 
r, and p gas ( r ) and Pth( r ) are the gas density and ther- 
mal pressure, respectively. These are measured within 
spherical shells, logarithmically-spaced in radius from the 
cluster center. 

We also measure the contribution of isotropic pressure 
from random gas motions, 



^rand = ^Pgas (°l + CT £ 2 



(2) 



where of and of are radial and tangential velocity dis- 
persions measured in spherical shells, using the peculiar 
velocity of the cluster dark matter within r$oo to define 
a rest frame. In the absence of radiative cooling, rota- 
tional velocities ar e negligible outside ~ 0.2r5oo, as was 



seen m 



Lau et al. (2009), however for completeness we 



treat their contribution to the effective pressure support 
separately. 



motions using approach similar to the treatment of the 
dynamical collisionless sys tem using the Jeans equation 
(Bin ney fc Tremaine|2 008). The total mass of the cluster 
within a radius r can be split into separate components 
corresponding to different terms of the Jeans equation, 



M tot (< r) = M HSE + M rand + M rot . 



(3) 



The streaming and cross terms are small, and we neglect 
them throughout this work. The rotational component, 
which contributes significantly to the pressure support in 
the core (< 0.5 rsoo), is given by 



M rot (<r) 



G ' 



(4) 



where v t is the mean tangential velocity. The component 
from random gas motions is given by 



rand 



(< r) = 



gas 



' d(p sag q?) 
dr 



G 



(5) 

Note that spherical symmetry for the gravitational po- 
tential and steady state are assumed in deriving the ex- 
pression for mass from the Jeans equation, and all the 
physical quantities at a given radius are averages over a 
radial shell. 

To avoid biases due to the p resence of substru cture and 
to maintain consistency with Lau et al.l (12009 ) , in mea- 



suring Mhse and P ran d we remove gas from within the 
tidal radius of subhalos of mass greater than 1O 12 /i _1 M . 
We further reduce small scale fluctuations in pressure and 
pressure gradient profiles by smoothing with a Savitzky- 
Golay filter. We have explicitly checked that neither sub- 
halo removal nor smoothing procedures have a significant 
effect on our results. 

It is important to note that all radial quantities in this 
paper, including Mhse, are measured within the corre- 
sponding "true" radius, ftme, as m easured directly from 
the simulations. |Lau et ah] ( |2009[ ) shows that if Mhse is 
instead measured within the radius as estimated from the 
HSE equation, r est , the hydrostatic mass bias withi n r^oo 
is increased by ~4% for relaxed systems (see also Nagai 
et al.|[2007b for results based on mock X-ray maps)^ 



3. EVOLUTION OF THE HYDROSTATIC MASS BIAS 

We begin by exploring in detail the evolution of a single 
object, CLIO, with a relatively quiescent merger history 
prior to a major merger at z ~ 0.72. This allows us to 
isolate the effect of a single merger and compare directly 
to previous studies tha t explo re idealized binary mergers 
(e.g., |Ricker fc Sarazin|2001| [Poole et al.||2006|) . 

The response ol the ICM thermal pressure during the 
merger is illustrated by Figure [l] Each frame is 2h~ 1 
Mpc comoving across and lh~ x Mpc comoving deep in 
projection, rsoo is marked with a thin white circle at each 
epoch, and r 2 oo — 1.5r 50 o and r 2 5oo - 0.5r 50 o. Initially 
CLIO is very relaxed, not having experienced a merger 
for 4.2 Gyr. Note the nearly circular pressure contours 
within rsoo in the first panel. The merging cluster, with 
mass ratio of 1:2, enters from the top of the frame with 
an orbit nearly in the plane of the projection, achiev- 
ing closest approach approximately ~ 0.5 Gyr after the 
start of the merger. This naturally divides the merger 
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Fig. 1. — Thermal pressure maps of CLIO at six epochs surrounding the last major merger. Each panel is 2h 1 Mpc (comoving) and 
l/i -1 Mpc (comoving) along the line of sight, centered on the most massive progenitor. The white circles depict the size of rsoo and the 
red dashed arcs mark the location of the merger shock. The time relative to the beginning of the merger (the epoch when the merging 
cluster cross r2oo) is shown in each frame. Frame (b) shows the epoch when the merging cluster crosses r^oo, frame (c) shows the time of 
closest approach, and frame (f) shows a time when CLIO has relaxed significantly after the merger. 
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Fig. 2. — The mass profile of CLIO at the same six epochs depicted in Figure [T] The hydrostatic mass estimate, Mhse, and the true 
mass, Mtme, are shown by the dashed and solid lines, respectively. The downward racing arrow marks the location of the incoming merging 
halo. The upward facing arrow marks the location of the shock wave created by the merger as it emanates outwards through the parent 
halo. 
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into three epochs. In the initial passage (panels a-c) 
the spherically averaged outer pressure profile becomes 
shallower due to the presence of the merging cluster and 
compression of the ICM between the two cluster cores. 
At the time of closest approach the inner pressure profile 
becomes strongly peaked, generating forward a nd reverse 



merger shock s, marked with a dashed red arc (Ricker 



Sarazin|2001 ). In this second phase of the merger (panels 



d and e), the merger shocks propagate outward, crossing 
rsoo at ~ 0.9 Gyr following the start of the merger. Af- 
ter several sound crossing times the merger shocks have 
passed out of the cluster central region and merge with 
the large-scale accretion shock. The last panel of Fig- 
ure [l] (f) shows CLIO at z = 0, - 7.2 Gyr after the 
merger. By this time the disruption of the merger has 
dissipated, leaving a relaxed, symmetric pressure profile 
similar to that seen at the start of the merger. 

Deviations from hydrostatic equilibrium can be eas- 
ily understood given the features seen in the pressure 
maps. In Figure |2j we compare the true mass profile 
(solid lines) to that predicted from the equation of hy- 
drostatic equilibrium (Equation [I] dashed lines) at the 
same six epochs depicted in Figure [I] The start and 
end states are again qualitatively similar, with a cluster 
core in near-equilibrium and progressively greater devia- 
tion with increasing radius beyond rsoo- The systematic 
underestimate of the total mass for relaxed clusters has 



been noted previously (|Rasia et al. 
|2007b[ [Q[u^ni1[2009t . 

The hydrostatic mass 



2004 Nagai et al. 



hydrostatic mass bias shows significant evolution 
in both radius and time. As the merging cluster makes 
its initial passage (panel b), with position denoted by the 
downward arrow, the true mass profile increases faster 
than the mass profile derived from the thermal pressure, 
temporarily increasing the (negative) hydrostatic mass 
bias at rsoo- 

The outwardly moving merger shocks, with approxi- 
mate positions shown by upward arrows in panels (d) 
and (e), produce locally very steep pressure gradients, 
and thus strong overestimates of the true mass. This is 
seen as a positive peak in the hydrostatic mass bias pro- 
file that propagates outward with the shocks. The width 
of these peaks is caused by a combination of spherical 
averaging and the smoothing procedure described in Sec- 
tion^ This evolution is in agreement with the evolution 
of the h ydrostatic mass profile during a merging event 
shown in |Puchwein fc Bartelmann| ( |2007| ) . 

The evolution ol the hydrostatic mass bias at a fixed 
radius is governed by the radial migration of these fea- 
tures across that radius. Broadly, the bias will decrease 
slightly during the initial phase of the merger. As the 
merger shocks move outward the bias will see a dramatic 
rise, peaking at shock crossing before falling once more 
to a large negative bias. This bias will diminish at all 
radii as the cluster continues to relax from the inside- 
out, however, at r > rsoo the bias persists to late times. 
As we will show in the next section, these features seen 
in the mass bias for the CLIO merger are generic to all 
the cluster mergers in our sample. 

We now expand our exploration of the hydrostatic 
mass bias beyond the time immediately surrounding the 
CLIO merger. The top left panel of Figure [3] shows the 
evolution in the mass bias at 3 different radii, r25oo? r 500? 
and r2oo, chosen for their relevance to X-ray and SZ ob- 



servations. The bottom panel shows the mass accretion 
history at r 5 oo normalized by the M 50 o at z = 0. A ver- 
tical line marks the beginning of the final major merger. 
The epochs shown in the previous figures are marked 
with red tick marks along the x-axes of the panels. 

The final major merger (depicted in the previous fig- 
ures) begins at ~6 Gyr. The mass of the cluster rises 
quickly as the incoming cluster merges with the parent 
cluster. About 1 Gyr after r2oo crossing, we see another 
sharp, very positive bias in the top panel. The bias falls 
again and remains negative through z = 0, oscillating 
slightly as the cluster relaxes. 

The evolution in the other radii is qualitatively simi- 
lar, also showing the short-lived, positive bias immedi- 
ately following the merger. Note that the peak at each 
radius occurs at a different time following the merger, 
with a delay of approximately ~ 0.5 Gyr between r25oo 
and rsoo and between rsoo and r2oo? reflecting the shock 
propagation we observed in panels (d) and (e) of Fig- 
ure [2] Immediately following the peak, the hydrostatic 
mass again underestimates the true mass and this nega- 
tive bias holds but steadily diminishes to z = 0. 

For comparison we also show in Figure [3] the evolution 
of two clusters selected from our full sample of 16. CL104 
is a massive cluster that forms from a complicated triple 
merger at z ~ 1.5 and remains mostly quiescent there- 
after. CL6 is intermediate in mass between CLIO and 
CL104, and experiences a nearly one-to-one merger at 
z « 0.6. Despite the differences in merger history and 
z = mass, the clusters show remarkable similarity in 
their evolution. All three clusters have peaks ~1 Gyr 
following the identified merger, and show similar relative 
offsets between the different radii. 

Cluster CL104 shows less prominent merger character- 
istics than seen in the other two clusters due to the triple 
merger, resulting in broader and smaller amplitude peaks 
at each radius. Following the major merger at z « 1.5, 
CL104 also undergoes a minor merger at z ~ 0.25. How- 
ever, this merger only has a mass ratio of 10:1. Con- 
sequently the minor merger does not impact the cluster 
enough to have a significant effect of the mass estimate. 

By z = 0, the mass bias for all three clusters ap- 
proaches zero. However, it is important to note that 
even at z = the hydrostatic mass bias within rsoo is 
still biased low by 5-10% in the outskirts. 

3.1. General Trends in the Hydrostatic Mass Bias 

The post-merger evolution in Figure [3] exists for all 
16 clusters in our simulated sample. This is shown ex- 
plicitly in Figure HJ where the average evolution in the 
hydrostatic mass bias is plotted for the entire sample. 
The final major merger is identified for each cluster (see 
Table [T]) and used to define a common point between 
clusters with varied accretion histories. The hydrostatic 
mass bias is measured at three separate radii, r25oo> ^500 ? 
and r2oo> linearly interpolated to equally spaced times 
relative to the last major merger, and averaged over the 
entire sample. Error bars denote the error on the mean 
for rsoo- The error bars for r25oo and r2oo are compa- 
rable, and omitted for clarity. Note that clusters that 
undergo low-redshift mergers will only contribute to the 
low-t 



merger 



portion of the figure. Consequently, while 
the behavior exhibited in the early epochs of Figure [4] is 
an average over the entire sample, the later epochs are 
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Fig. 3. — Top: The time evolution of the mass bias (top) for three clusters, CLIO, CL104, and CL6 at three radii r25oo (red, dotted), rsoo 
(black, solid), and r2oo (blue, dashed). Bottom: The mass accretion history of the three clusters at rsoo, normalized by M500 at 2 = 0. 
The vertical gray line marks the beginning of the latest major merger for each of the clusters. The red ticks mark the epochs in CLIO 
corresponding to the panels in Figures 1 and 2. 

necessarily biased to a small number of early forming ob- 
jects. We plot the evolution of the mass bias to 9.25 Gyr 
following the last merger at which point fewer than three 
clusters contribute to the average. 

The qualitative features for the individual clusters pre- 
sented in Figure [3] are also apparent in the sample- 
averaged evolution in Figure [4] Immediately preceding 
the merger, the hydrostatic mass underestimates the true 
mass by ~ 10% at all radii. This bias grows to nearly 
30% in the outskirts shortly after the merger begins. Fol- 
lowing the peaks associated with the merger shocks, the 
hydrostatic mass bias decreases over the subsequent 8 
Gyr from -15-20% to 3%, 10% and 18% within r 250 o, 
rsoo and r2oo, respectively. 

The peaks caused by the propagating merger shocks 
have smaller amplitude and are broader than the indi- 
vidual cluster histories seen in Figure [3j Therefore, the 




mean mass bias remains negative at alTt merger . This is 
caused by a variety of factors, including differences in 
the merger mass ratio and the orbital impact parameter, 
which changes the time offset between the peak and the 
start of the merger. Since the peaks are narrow in time 
(typically less than 0.5 Gyr), and we are primarily con- 
cerned with the average evolution after the merger, we 
make no attempt to reduce this broadening. 

We also examined how the mass bias depended on var- 
ious parameters, such as merger ratio, final cluster mass 
and redshift. We found little dependence on present mass 
of the cluster, with the larger clusters experiencing a 
slightly greater peak in the mass bias immediately fol- 
lowing merger. 



4. 



INFLUENCE OF NON-THERMAL PRESSURE ON THE 
HYDROSTATIC MASS BIAS 



Figures [3] and [4] demonstrate that the magnitude of 
the mass bias decreases steadily following a merger, yet 
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Time since last merger (Gyr) 



Fig. 4. — Averaged mass bias as a function of time elapsed since 
last merger (in Gyrs) for the sixteen clu sters. A more detailed 
discussion of this figure is located in §3.1| The biases are plotted 
at radii r25oo (red, dotted), rsoo (black, solid), and r2oo (blue, 
dashed). The error bars show the la error on the mean at rsoo- 



never f ully disappears. If, as has been previously sug- 
geste d ( |Lau et al.||2009| |Nagai et al.||2007b| jRasia et al. 
2004), residual motions in the 1C1V1 account tor this bias, 
the non-thermal component of pressure due to these mo- 
tions should see a similar evolution to the hydrostatic 
mass bias. We now investigate the evolution of this non- 
thermal component in more detail. 

The top panels of Figure [5] show the evolving contribu- 
tion of random motions to the total ICM pressure (left 
panel) and pressure gradient (right panel) averaged over 
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Fig. 5. — Top: Fractional contribution of random motions to the total effective pressure (left) and pressure gradient (right). Bottom 
left: The components of the total ICM energy. Bottom right: Averaged corrected mass bias as a function of time elapsed since last merger 
(in Gyrs) for the sixteen clusters. In all figures, the quantities are plotted at radii r25oo (red, dotted), rsoo (black, solid), and r2oo (blue, 
dashed). Our sample is divided into relaxed and unrelaxed clusters at ^merger = 4 Gyr (marked by the vertical dashed line). The error 
bars show the la error on the mean at rsoo- 



all 16 clusters in the same manner as in Figure [4] The 
measurement of t he no n-thermal pressure support is de- 
scribed in Section 2.3 During the merger, the effective 
pressure from random motions is ~ 30% at rsoo and has a 
strong trend with radius, increasing from ~ 22% at r25oo 
to ~ 38% at r2oo- This radial trend persists through- 
out the merger and during the following relaxation. The 
non-thermal pressure fraction drops rapidly over the next 
1-2 Gyr, and steadily decreases thereafter from ~ 20% to 
~ 10% at r 50 o. 

The radial pressure gradient shown in the top right 
panel of Fig. [5[ an analog to the thermal pressure gradi- 
ent in the equation of hydrostatic equilibrium (Eq. [I]), 
shows similar trends in both radius and time since 
merger. Initially, random velocities contribute ~ 20%, 
~ 26%, and ~ 32% to the non-thermal pressure gradi- 
ent for r25oo, ^500, and r2oo, respectively. This falls to 
between 2 — 12% after 8 Gyr, in line with the residual 
hydrostatic mass bias at those radii. 

Our findings are consistent with the measurements of 
previous work. |Lau et ah] ( |2009[ ) also sees small scatter 



in the gradient of non-thermal pressure at rsoo from 4- 
10% among clusters identified as relaxed at z = and 
a larger range of 3-18% for unrelax ed clusters (see their 
Figure 3). While [Lau et al.| ( |2009| ) used the same sim- 
ulations as this work, our mass biases are also in gen- 
eral agreement with the me asurements from o t her simu- 
lations, e.g. GAD GET- 2 in|Rasia eFaT ](|20 p6[);|Piffaretti 
fc Valdarnini| ( |2008j ); [Battaglia et al.| ( |2Ull[ ). The strong 
radial dependence m the non-thermal pressure fraction, 
with larger non-thermal s upport at greater radii is con - 
sistent with the results by |Iapichino fc Niemeyer| (2008), 
|Battaglia et al.| ( |2Q11| ), and |Vazza et al.| ( |2Ull| ). In ad- 
dition, the dependence on dynamical state seen in top 
panels of Fig ure [5] is consistent with the findings of Vazza 
|et al.| ( |2Q11| ). They also demonstrate that clusters un- 
dergoing mergers have a significantly larger fraction of 
non-thermal pressure support than relaxed or relaxing 
systems. 

We expect the decay of the non-thermal component 
seen in the upper panels of Figure [5] is due to the conver- 
sion of random velocities to thermal energy via shocks 
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and numerical dissipation (no explicit physical viscosity 
was included in our simulations). Demonstrating this 
explicitly is difficult, however, as the total energy in the 
ICM increases by continued accretion subsequent to the 
merger, as can be seen in the bottom panels of Fig- 
ure [3J To account for this growth we rescale the energy by 
the self-similar mass and redshift scaling M 5 / 3 £ , (z) 2 / 3 , 
where E(z) — H(z) /Hq is the dimensionless Hubble con- 



stant (Kaiser 1986[_Nagai et al. 2007a). In the bottom 



left panel of Figure jo] we show the the evolution of the to- 
tal (solid) , thermal (dotted) and random kinetic (dashed) 
energy of the ICM within rsoo, scaled by the self-similar 
expectation. To show the evolution of energy at differ- 
ent radii, we also plot the random kinetic energy within 
three radii: r25oo (red), rsoo (black), and r2oo (blue). 

The total energy content of the ICM remains essen- 
tially constant beyond the 2 Gyr period following the 
merger. The time surrounding the merger shows a sig- 
nificant decrease in E tot , however this is due not to a 
decrease in energy but rather the temporary increase in 
M500 caused by the presence of the merging cluster in 
the cluster core. Much of this mass will move to larger 
radius, decreasing the mean density and thus M500, as 
can be seen in the accretion histories shown in Figure [3] 
The apparent decrease in energy is a reflection of this 
deviation from self-similarity, rather than a change in 
the total ICM energy. Since the total energy remains 
constant, we can associate directly the decrease in ran- 
dom kinetic energy with the increase in thermal energy. 
This therma l izatio n is in agreement with the findings of 



Vazza et al. (2011) who show that the fraction of turbu- 
lent energy to thermal energy increases with decreasing 
time since major merger. This evolution in the thermal 
energy has important implications for the behavior of 
cluster mass proxies such as the integrated SZ flux (e.g. 
|Nagai||2006| IShaw et al.|2008l |Yang et al. 2010) or X-ray 
temp erature (e.g Ritchie & Thomas 2002; Rowley et al 
2004; |Rasia et al.||201l| ) which" are sensitive to the ther 
mal energy content of the ICM, and will be explored in 
greater detail in future work. 

It is interesting to ask at what point the large bulk 
velocities associated with the merger progenitors reach 
equilibrium with the new gravitational potential. We 
produce a corrected mass estimate, M tot , using Equa- 
tion [3] and measuring the terms due to random motions 
and rotation. 

The lower right panel of Figure [5] shows the evolution of 
M to t following the merger. The period up to « 2 Gyr fol- 
lowing the merger shows rapid evolution, with the mass 
bias ranging from a positive 20% bias to a negative 10% 
bias. This is driven primarily by the evolution in Mhse 
caused by the outwardly propagating merger shock seen 
in Figure [I] The mean M to t after 4 Gyr is essentially 
unbiased, with a mean of 0.3% and « 8% scatter at 
rsoo- This indicates that the velocities at this radius 
have reached equilibrium with the cluster potential. 

4.1. Correcting the Hydrostatic Mass Bias at z = 

The dependence of M tot on the time since major 
merger suggests a promising division between "relaxed" 
and "unrelaxed" clusters may be obtained by splitting a 
sample between objects that have and have not experi- 
enced a major merger within the past 4 Gyr. This results 
in a sample of 10 relaxed clusters at z — (see Table [lj. 
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Fig. 6. — The averaged mass measurements for the relaxed sub- 
sample, split at tmerger of 4 Gyr at z = 0. We show Mtot (solid 
black) and the mass corrected using the relation from [Shaw et al~] 
(2010) fit to our subsample, M CO rr (dashed red). See text tor de- 
scription of these mass calculations. Mhse is shown with the solid 
red line for comparison. 

Note that this includes all 6 clusters in Nagai et al. 07 
identified as relaxed clusters based on morphological clas- 
sification. Our relaxed subsample also contains an addi- 
tional 4 clusters that were missed by the morphological 
classification scheme. This suggests that these clusters, 
while dynamically relaxed, appear morphologically dis- 
turbed in the X-ray images. 

Figure [6] shows the average Mhse (solid red line) and 
Mtot (solid black line) radial profiles for our relaxed clus- 
ters at z = 0. The inclusion of M rot and M ran d in the 
calculation of M to t removes the bias that is apparent for 
Mhse out to « 1.5r 5 oo, with a slight increase in scatter 
from 8% to 11% at r 50 o and 9% to 10% at r 20 o. 

Since P ran d is not accessible observationally, we pro- 
pose a model based on the average P ran d profile mea- 
sured in our simulations. We neglect the contribution 
from rotation, which is < 2% of the mass correction 
beyond 0.7r5oo, and also assume velocity anisotropy, 
P = 1 — cr^/2cr^, is zero. While j3 is nonzero and increases 
radially, we find that the inclusion of (3 as measured from 
the simulations has little effect on the accuracy of our 
model. We parameterize the fractional contribu tion of 
the non-thermal pressure as in Shaw et al. (2010), 



rand 



rand 



th 



a(r/r 50 o) r 



(6) 



Using the average P ran d and P tn profiles for our relaxed 
sample and fitting over the radial range 0.2 < r/r 5 oo < 
2.0, we obtain best fit values a = 0.14 ± 0.09 and 
n = 1.16 ± 0.17. Results are mildly dependent on red- 
shift in the re dshift range z < 0.2 (see |Shaw ctllpQIQ 



Battaglia et al. 2011 for redshift evolution). These dif- 
fer from the best fit values given in Shaw et al. (2010) 
of a = 0.18 and n = 0.8, which were derived for the 
complete sample, including unrelaxed clusters. The re- 
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TABLE 2 

Mean Biases in Mass Estimates at z=0 



Mass Estimate 


Sample 1 


Bias ± Scatter 2 


^500 


^200 


M HS e 


Full 


-0.10 ±0.16 


-0.06 ±0.19 


M HS e 


Relaxed 


-0.07 ±0.08 


-0.14 ±0.09 


Mtot 


Relaxed 


-0.01 ±0.11 


-0.01 ±0.10 


Af C orr 


Relaxed 


0.01 ±0.09 


0.00 ±0.11 



1 Relaxed clusters are clusters that have not experienced a major 
merger in at least 4 Gyr. 

2 The error on the mean is given by the scatter divided by — 1, 
where N is a number of clusters. 

suiting corrected mass profile contains only observable 
quantities p gas , P t h, and dP^/dr, and is given by 



1 



dR 



th 



Gp gas \1- a(r/r 500 ) n dr 

an(r/r 50 o) n_1 



nsoo(l - a(r/r 50 o) n ) 2 



th 



(7) 



The dashed red line in Figure [6] shows the average 
M corr for the sample of relaxed clusters. Using the av- 
erage non-thermal pressure profile we are able to re- 
cover the true mass between 0.7-2r5oo- At smaller radii, 
r < 0.7r 5 oo, the neglected rotational component becomes 
important and the mass is again underestimated. A com- 
parison of the mean mass biases for M corr and Mhse at 
two radii are shown in Table [2] It is worth noting that the 
mean mass biases for both the full and relaxed subsampl e 



are in good agreement with those of Lau et al.| (120091) 
despite differences in gas physics employed in the respec- 
tive runs; i.e., non-radiative simulations in this work and 



the simulati ons including cool and star formation in Lau 
et al. (2009). The M corr mass estimate is able to recover 
the true mass on average to within 1% with 10% scatter. 
This is a significant improvement over the mean HSE 
mass bias as measured for the full sample of 5-10% with 
15-20% scatter. 

It is important to note that this model is not appli- 
cable to unrelaxed clusters, for which the fundamental 
assumptions in the hydrostatic equilibrium equation do 
not hold. However, for the full sample, we find that 
the mass bias is nearly constant with very little radial 
depe ndence, which is in a greement with previous work 
(e.g. |Battaglia et al.||2011| ). Applying a constant mass 
correction ol 107o (equivalent to using our model with 
parameter valu es a = 0.1 and n = 0), as has been done 
previously (e.g. Mantz et al. 2008), will remove the bias 
on average for a sample including both relaxed and un- 
relaxed clusters. The scatter in this correction is large, 
however, at 15-20%. Therefore, our correction, though 
only applicable to a subsample, can provide significant 
improvements to the cluster mass measurements when 
such a subsa mple is applicable, such as gas fraction mea- 
surements QAllen et aL||2008| [LaRoque et al.||2006| ). We 
note that the correction should in principle be a function 
of redshift, as the contribution of non-thermal pressure 
grows to high er redshift when clusters are dynamicall y 



younger (e.g. Battaglia et al. 2011 |Shaw et al. 2010). 



However, since we do not sample the lull accretion histo- 
ries of clusters at all redshifts, we report best fit values 
for our model only at z = 0. 



5. DISCUSSION AND CONCLUSIONS 

Using high-resolution adaptive mesh refinement cos- 
mological simulations of a sample of sixteen clusters, we 
have performed an analysis of the effect of major mergers 
on the hydrostatic mass estimate. Rather than examin- 
ing the properties of galaxy clusters at a single epoch, we 
utilize merger trees to follow the evolution of individual 
objects during and after mergers. We focused specifically 
on the characteristic features of the hydrostatic mass bias 
during the merger and the subsequent evolution of the 
random ICM motions produced. 

Our principal results can be summarized as follows: 

• Major mergers are the dominant sources of bulk 
motions providing non-thermal pressure support in 
the ICM and therefore bias in the hydrostatic mass 
estimate. 

• Well after a major merger, the hydrostatic mass 
estimate underestimates the true mass within rsoo 
by 6-10%. 

• During a major merger, random kinetic motions 
contribute ~30% of the total ICM pressure at r 5 oo- 
Subsequently, the energy in these motions are con- 
verted to thermal energy over time, comprising a 
residual 10% contribution when the cluster is well 
relaxed. 

• We provide a method for correcting the hydrostatic 
mass estimate for relaxed clusters based on the 
thermal pressure profile that recovers the true mass 
within ~10% from 0.7-2r5oo- 

In detail, we find the evolution of the hydrostatic mass 
bias can be divided into two epochs. During the merger, 
the infall of the smaller structure generates asymmetrical 
shocks that emanate from the center of the parent clus- 
ter. These shocks result in a mass bias peak that prop- 
agates outward over the course of a Gyr until it reaches 
the edge of the cluster. 

As the clusters relax, the hydrostatic mass bias evolves 
more slowly, decreasing from ~20% to 6-10% within rsoo- 
Even after 8 Gyr, the mass bias in the innermost regions 
(< ^2500 ) fluctuates around 3%. The mass bias increases 
with increasing radius, reaching 18-20% at r2o o, which is 
in ag r eement with previous simulations (e.g. , |Kay et al. 
2004; |Rasia et al.] [20061 INagai et afl|2007b| [Lau et al 
2009b and o bservations (e.g., |Mahdavi et al.|2008[ |Zhang 
|eTaL|2010| ). 

The evolution of the hydrostatic mass bias in the sec- 
ond epoch, after the merger, is driven by non-thermal 
pressure support from bulk motions. The random ki- 
netic energy decays after the merger, thermalizing into 
the surrounding ICM. Despited continued accretion, the 
fraction of total ICM energy in residual bulk motions falls 
from ~10% immediately following the merger to ~5% 
long after. Partly this reflects the inability of smoothly 
accreted material to excite new motions in the ICM and 
partly the thermalization of existing random motions. 
The increase in the thermal energy drives the observed 
decrease in the hydrostatic mass bias. 

Since the bulk motions in the ICM drive the hydro- 
static mass bias, we investigated methods for improv- 
ing the mass estimate by accounting for the non-thermal 
pressure provided by random gas motions. 
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Within 4 Gyr of a major merger, a cluster has not yet 
virialized and is not in equilibrium. However, after this 
time, a cluster has sufficiently relaxed that with accurate 
treatment of the random pressure, we are able to recover 
the true mass. We therefore use this criteria to split 
our sample into "relaxed" and "unrelaxed" by dividing 
into clusters that have and have not experienced a major 
merger in the last 4 Gyr. 

We provide a method for correcting the hydrostatic 
mass estimate for relaxed clusters by harnessing the sim- 
ilarity in the evolution of the non-thermal pressure pro- 
files. This method combines the true thermal pressure 
profile for each cluster individually with the average non- 
thermal pressure profile of the relaxed sample. At z = 0, 
our model recovers the true mass to within ~10% in the 
radial range 0.7-2 rsoo- 

There are a few issues to bear in mind when interpret- 
ing our results. First, our simulations neglect radiative 
cooling, energy feedback from stars and active galactic 
nuclei, and additional sources of non-thermal pressure 
such as cosmic-rays and magnetic field. While these ef- 
fects are critical for cluster cores, their effects on cluster 
outskirts remain uncertain and must be explored in the 
future studies. Uncertainty in our results may also arise 
from the lack of viscosity and our treatment of the ICM 
as a single fluid, neglecting various effects important for 
plasmas. Finally, our measurements use full three dimen- 



sional information from the simulations averaged using 
spherical symmetry and therefore do not account for pro- 
jection effects. We exclude gas contained within subhalos 
of mass greater than 10 12 h~ x M , which will differ some- 
what from observationally identified peaks in emission. 
These differences may account for our somewhat smaller 
mass bia ses than those measured from mock X -ray anal- 
yses (e.g. Nagai et al.|2007b||Rasia et al.|2004[ ). However, 
despite these caveats, we are confident that it is possi- 
ble to gain better control on the bias in the hydrostatic 
mass estimate of clusters by taking into account its de- 
pendence on dynamical state. In the future work, we will 
explore possible morphological classification schemes for 
defining these subsamples based on observational data. 
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